Notes on the Qubit Phase Space and Discrete Symplectic Structures 
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We start from Wootter's construction of discrete phase spaces and Wigner functions for qubits 
and more generally for finite dimensional Hilbert spaces. We look at this framework from a non- 
commutative space perspective and we focus on the Moyal product and the differential calculus on 
these discrete phase spaces. In particular, the qubit phase space provides the simplest example of a 
four-point non-commutative phase space. We give an explicit expression of the Moyal bracket as a 
^ ' differential operator. We then compare the quantum dynamics encoded by the Moyal bracket to the 

^ , classical dynamics : we show that the classical Poisson bracket does not satisfy the Jacobi identity 

thus leaving the Moyal bracket as the only consistent symplectic structure. We finally generalizes 
' our analysis to Hilbert spaces of prime dimensions d and their associated d x d phase spaces. 
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It is well-knovifn that standard quantum mechanics can be mathematically reformulated in almost classical terms 
] using Wigner functions and the Moyal product. This is achieved through the Weyl transform, which maps phase- 
space functions to Hilbert-space operators and vice-versa. A by-product is a reformulation of quantum mechanics 
as a non-commutative geometry with a non-commutative (Moyal) product between functions on the phase space. 
This construction was initially introduced for systems evolving in a continuum and thus for infinite dimensional 
' Hilbert spaces such as L^(R"). More recently, Wootters and collaborators introduced a similar construction for finite 
dimensional Hilbert spaces with a discrete phase space and discrete Wigner quasi-probability distributions [E 0, Hi 
(see also Q). For instance, the qubit with its two-state Hilbert space is described in term of a four-point phase space 
^ [ presented as a 2 x 2 lattice. In the present notes, we do not introduce a new phase space construction or discuss the 
^ r use of Wootters' formalism in quantum tomography or related applications. We are interested in Wootters' framework 
' from a point of view of non-commutative geometry and more specifically we focus on the properties of the discrete 
T-H , Moyal product and the related symplectic structure. Indeed this construction provides very simple examples of non- 
' commutative spaces with finite number of points. These are provided with discrete non-commutative ★-products, 
[ which reflect the non-commutative product between operators on the Hilbert space. 
\^ • Starting with the single qubit, we introduce its 2 x 2 lattice phase space and its discrete Moyal product. We define a 
I differential calculus on his discrete space and provide an explicit formula of the Moyal bracket in term of the discrete 
■ differential operators. This shows that this phase space is equipped with a discrete symplectic structure. We further 
compare this Moyal bracket with the natural Poisson bracket defined in term of the introduced differential calculus. 
They both have a very similar expression but they lead to substantially different dynamics: the Moyal bracket leads 
to negative probabilities as expected since it describes quantum dynamics, while the Poisson bracket always keep 
^ . probabilities positive. However, this Poisson bracket is shown not to satisfy the Jacobi identity (associativity) and 
' thus does not provide a consistent symplectic structure. Somehow, this is not surprising since the qubit is a true 
quantum system with no classical analog: there is no ?i — > limit in which the Moyal product would be approximated 
by a classical Poisson bracket on the discrete phase space. 

In a second step, we generalize our analysis to Hilbert spaces with higher dimension. Focusing on the case of a 
prime dimension d, we follow Wootters' construction and we define the discrete Weyl map and the Moyal product on 
the discrete phase space. We give explicit formulas for the Moyal product and brackets. We compare these expressions 
with the standard Poisson and Moyal brackets on the usual two-dimensional continuous {q, p) phase space (for a one- 
dimensional system). We show that the continuous Moyal bracket evaluated on certain holonomic observables gives 
exactly the discrete Moyal bracket: the discrete phase space is thus an exact truncation of the standard phase space 
where the space of functions over the phase space is restricted to (the Hilbert space generated by) these holonomic 
observables. This correspondence also shows that the Planck constant h somewhat goes effectively as the inverse of the 
dimension d~^. Finally, we define the discrete differential calculus, discuss the free evolution with a "p^"-Hamiltonian 
and identify canonical transformations leaving the symplectic structure invariant. Once again, we conclude that the 
discrete Moyal bracket defines the only consistent symplectic structure and that there does not seem to be any analog 
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of a classical Poisson bracket expressed in term of the discrete differential calculus and satisfying the Jacobi identity. 



I. A QUICK REVIEW OF THE QUBIT PHASE SPACE 

We work with a qubit, i.e living in a spin-^ representation of SU(2). The Hilbert space is two-dimensional. We 

call X,Y^Z the three Pauli matrices ^, satisfying the usual Lie algebra eommutator [^,1^] = 2iZ. These operators 
X, Y, Z have eigenvalues ±1. We denote as usual | T); I i) the two eigenvectors of Z, and we introduce the eigenvectors 
I ), I — >). The projectors on these eigenvectors read: 

|T)(TI = ^(/ + ^), = + 

\i){i\ = \{i-z), \^){^\ = \{I-X). (1) 

We define the following four (Hermitian) observables: 

D_+ = \{I + X-Y-Z), D++ = \{I + X + Y + Z), 

D__ = ^{I-X + Y-Z), D+_ = ^(I-X-Y + Z), (2) 

or more compactly D^^i = {I + e'X + ee'Y + eZ)/A. These -D-operators form a POVM and define an abstract four 
point space ±±. This defines the discrete phase space for the qubit. This construction has the property that the 
sums of the iJ-operators along the lines of this space give the projectors on the eigenvectors of X and Z: 

D__ + D_+ = \l){l\, D__+D+_ = |H(^|, 
2?+_+D++ = |t)(T|, D_+ + D++ = \^){^\. 

We also have _D__ + -D++ = (/ + Y)/2 and D |- + = (/ — Y)/2 along the diagonal lines. In that sense, the 

operator represent the point (|, <— ) and so on. Since we can not diagonalize simultaneously X and Z, we can not 

fully distinguish -separate- these four points. This translates to the fact that the D operators do not commute with 
each other, e.g {D , = i{X + Y)/8. 

This four point space is our discrete phase space. The eigenvectors J,, t define the "q" coordinate direction while 
■*— , define the "p" momentum direction. 

These D operators satisfy the following trace identities: 

Va,/3=(±±), tr£>„ = ^, tr£)„-D/3 = ^^a/?- (3) 

Resulting, the D's form an orthogonal basis for linear operators on H and we can decompose any operator in this 
basis. We have : 

F = 2 ^ f„D„, with /„ = tr FD^, (4) 

a 

for all operators F on H. The factor 2 is simply the dimension of the Hilbert space H. Notice that since J2a -^a ~ ^' 
we have trF = J2a /« functions on the phase space. 



^ Our convention for the Pauli matrices are: 



--{u)^ --(;-°,) 
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We can write density matrices in this basis. A density matrix p, satisfying = p,tv p = 1, can be naturally 
decomposed on the Pauli matrices: 

p=^{I + aX + bY + cZ), a,b,ceR. 

Since det p= (1 — — 6^ — c^c)/4, the positivity requirement p > means that + 6^ + < 1. We then define the 
projections Pa = tr DaP : 

1 1 
P-+ = -{1+a-b-c), P++ = -{l + a + b + c), 

p— = ^{I-a + b-c), p+_ = ^{1-a-b + c). (5) 

These are the probabilities that the system is located on each of the four points (±±) of our discrete phase space. 
We easily check that these probabilities sum up to 1. 

Finally, we can define a notion of "classical" states p for which all our projections are positive, Va, pa > 0. 



II. STAR PRODUCT, DIFFERENTIAL CALCULUS AND SYMPLECTIC STRUCTURE 

Wc define the multiplication between functions (obscrvablcs) on the discrete phase space as the operator multipli- 
cation between operators on H. More precisely, for functions fa = tiDaF and ga = tr DaG, we define: 

if^g)a = ltvDaFG. (6) 

The factor i is simply due to the factor 2 entering the map between operators and functions on the phase space. It 
ensures that the constant function = 1 is the unity for the ★-product. This product is obviously associative but 
not commutative and we define our discrete Moyal bracket as its commutator : 

{{f,gh)a = ~tTDa[F,G]. (7) 

The fi-factor is introduced by hand so as to match the expression of the symplectic structure in term of the Moyal 
★-product in standard quantum mechanics ^. This bracket is obviously anti-symmetric and satisfies the Leibniz rule 
and the Jacobi identity by construction (since the operator commutator does) ^. 

We give the simple examples of the operators X and Z. The projections of Z on the Da's are : 

Z__=tr£>__Z = -^ = Z_+, Z+- = Z++ = +^. 

Therefore the operator 2Z defines the coordinate g = ±1 in the phase space. Similarly, we have X = X-| = —1/2 

and = = -1-1/2, so that the operator 2X defines the momentum coordinate p = ±1. We can compute their 
Moyal bracket: 

{{q,ph)a = -jtTDa[Z,X] = ^trDaY. 



^ We remind the reader that, in standard quantum mechanics, Moyal bracket only matches the Poisson bracket at first order in h. 
^ For the matrix commutator, the Jacobi identity reads simply: 



[F, [G, H]] + [G, [H, F]] + [H, [F, G]] = 0. 
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The trace of D^Y gives ^ at the points a = ( ), (++) and — i on the points ( — h), (H — ). Therefore, we obtain the 

foUowing bracket: 

{q,ph = 2|. (8) 

Let us stress that the Planck constant H enters directly in the Moyal bracket because we have put it by hand in our 
definition ([7]) above. More interesting, we point out that such a quadratic structure for a Poisson bracket is usually 
the signature for a quantum group (and it gives directly its i?-niatrix). 

We now introduce differential operators. Following a standard strategy in non-commutative geometry, we define 
the derivation with respect to an operator 0, or equivalently with respect to the corresponding function 0, as the 
commutator with O : 

idefh^-ttTD4e,F], or equivalently def^2h{9JU. (9) 

Oq obviously satisfies the Leibniz rule d{f -kg) — f -k dg + df * g since the commutator [•, •] does. For instance, we 
define the derivative 8/ dp along the p-direction by considering the operator = Z. A short calculation yields: 

(dzf)— = (/++ - /+-), idzf)-+ = -(/++ - /+-), 

{dzfh- = -{f-+-f-), {dzf)++ = {f-+-f-)- (10) 

Similarly, dx defines the derivative along the g-direction and the derivation with respect to the K-operator gives 
variations along the diagonal lines Looking at the expression of the derivative dz, we notice that these formulas 

are not exactly as expected since (dzf) computes the finite variation in the p-direction at the point (H — ) instead 

of at the point ( ). To understand this point shift, we look more closely at the translations defined on our discrete 

phase space. 

For instance, Z generates translations along the p-direction : 

ZD^^Z = D^+, ZD^+Z = D^^, ZD+^Z = D++, ZD^^Z = D+^, 

where we remind that Z^^ = Z. Thus we can define an operator Tz acting on a function / as {Tzf)a = tr DaZFZ = 

tr ZDaZF. This gives (Tzf)-- — f \- and so on. This allows to define another differential operator computing 

functional variations in the direction generated by Z naturally as Szf = {Tzf — /): 

{Szf)a = trDaZFZ - ivDaF = tr Da[Z,F]Z = - tr Da[Z, FZ]. (11) 

The relation between the two differential operators dz and Sz is simple: dzf is equal to 6zf, with / corresponding 
to the operator FZ while / is given by F. This explains the point shift noticed early which is due to that extra- 
multiplication by Z. The new operator Sz does not satisfy the same Leibniz rule as the differential operators Oq, but 
we obtain a (slightly) deformed Leibniz rule: 

Szifg) - {Tzf)kSzg + Szfkg. 
In the following, we will denote the difference operators in the q and p directions respectively as 5q = 5x and Sp = Sz- 

We are now ready to give a full explicit expression for the ★-product: 

(/ * 5)a = ^ tr D^FG = 2 ^ ff,g^ tr{D^DpD^). 

/3,7 



We compute 

=(/-+-/++), (9x/)--+ = -(/— -/+-), (9x/)+- =-(,f-+ -/++), (9x/)++ = (/—-/+-), 
(c>y/)__ = (/_+-/+_), = (/++-/_-), =-(/++-/__), (9y /)++ = -(/_+ -/+_). 
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We can compute explicitly the traces of cubic polynomials of the f-matrices. For a = (++), we give 16tr(DQ,Z)^Z)^) 
in the following table with /3 labeling the columns and 7 the lines: 





++ 


+- -+ 




++ 


5 


1 1 


1 


+- 


1 


1 (-l-2i) 


(-l + 2i) 


-+ 


1 


+ 1 


(-1 - 2i) 




1 


(-l-2i) (-l + 2i) 
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This leads to the following formula: 

5 

if* 9)++ = g/++5++ 
1 



+ if++9+- + f++9-+ + f++9— + f+-9++ + f-+9++ + f—9++ + f+-9+- + f-+9-+ + f—9—) 



1 



1 



+- (-1 + 2i) {f+-9-+ + f-+g— + f—g+-) + -(-1 - 2i) {f+-g— + f—g-+ + f-+g+-) • 



(12) 



We first notice that even if the original functions fa and ga are real, their product has a priori a non-trivial imaginary 
part. This product can actually be written as /++(/++ + [difference terms]. The symmetric terms are real while the 
anti-symmetric terms are purely imaginary and contribute to the Moyal bracket. More precisely, we have: 



{{f,9}*)++ = —{f+-g-+ + f-+g— + f—g+--f+-g—-f—g-+-f-+g+-), 



1 
1 

2h 



(/+_ - f-){g-+ - g-) - (/-+ - f-){g+- 

Sqf^p9 - ^pf^q9 



(13) 
(14) 



We notice that we recover the dqAdp structure of the usual Poisson bracket, except for the point shift (++) ( — ). 
In fact, one can show ^: 



Va, tr D„ [F, G] = i tr DaY{XFX - F)Y tr DaY{ZGZ - G)Y - {F ^ G) 
which translates to the more compact formula in terms of functions on the (discrete) phase space^: 

{/> 9}* = -^'^Y [SqfSpg - SpfSqg] , 



(15) 



(16) 



where the product between Sf and 5g is the standard (point) product and not the ★ product. As we will see below, 
the Ty translation is very important when considering the dynamics on this phase space : it seems to be the difference 
between our Moyal ^-bracket for the qubit phase space and its classical counterpart. We could not find such a simple 
formula in term of the other differential operators dp^q. 



Remember that Ty = TzTx- 

We first prove the formula ti DaZ = 2tr DaX ti DaY which holds for all permutations of X,Y,Z and which follows from the simple 
calculations: 

Va = (ee')> 2trD,,/X = e', 2tr D^^,Y = ee' , 2tr D^^iZ = e. 
Then it is straightforward to show that: 

Va, tr D„ [F, G] = i [ tr {ZFZ - YFY) tr Da {XGX - YGY) - (F <-» G) ] , 

by checking it on the operator basis F. G G {I , X,Y, Z}. More precisely, this identity is trivial satisfied when F oi G arc / and when 
F = G. Moreover, since the formula is anti-symmetric, we only need to check explicitly three cases, {F,G) = {X,Y), [Y, Z),{Z, X). 
Finally this equation is equivalent to the expression given above. 
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III. EVOLUTION: QUANTUM VS "CLASSICAL" 



A Hamiltonian evolution dp/dt — —jr [H, p] can be translated in terms of functions on the phase space and ★-product: 

f = {H,pK. (17) 

For a single qubit, the only non-trivial (Hermitian) Hamiltonian (up to unitary transformation and multiplication by 
a constant factor) is H — X. As we showed earlier, this corresponds to the function p on the phase space and thus 
should simply generate motion along the q-axis. Using the previous calculation of dx, we get the following equations 
of motion: 

^ = (p-+-p++), ^=.(p++-p_+). (18) 

We check the conservation of probabilities along the classical trajectories: 

d{p++ + P-+) _ Q _ d(p+_+p__) 



dt dt 

where a trajectory is given by a fixed p (the second sign) and the two different signs of q. This translates the fact 
that eigenstates of X are stable under the evolution. 

We can check that these equations of motion naturally lead to negative probabilities. Physically, this is a signature 
of a quantum regime. Mathematically, this is due to the point shift in the discrete derivative dx- 

To solve the system, we call A the constant of motion: p++ + p |_ = A and p+- + p = 1 — ^1 (since trp = 1 is 

fixed). Inserting this in the equations, we obtain a second order equation for, say, p++: 

^-f±t = A + l-Ap++, (19) 

which finally get solved by p++ ~ K cos{2t) -f (A + l)/4 where K is another constant of integration (there could be 
an oscillation in sin too). Similarly, we derive: 

3A — 1 A — 1 3 — 3A 
p_+ = -A'cos(2t) H — , p = -Ksm{2t)^ — , p+^ = K sm{2t) ~\ — . 

We check that these probabilities sum to 1 as wanted. Then, due to the oscillatory behavior, we can easily get negative 
probabilities. For instance, taking A—1 and K — 1/2, we have p+^ — 1 and all other probabilities equal to at the 
initial time t — 0, but p becomes negative as soon as the system starts to evolve. 

Now, we can compare this dynamics with the "classical" evolution defined by the following equations of motion: 

dp++ , . dp^+ 

= (P-+ - P++), —j^^{P++-P-+), 

dp+- , . dp-- , . 

^ = (P+--P-), -^^{p-~P^-). 

Such an evolution would keep probabilities positive all the time (if they are positive to begin with). That's why we 
call "classical" this choice of dynamics. The main difference between these classical equation of motion is that the 

two trajectories (++ ^ — h) and ( — V ^ ) are completely decoupled, while the true quantum motion actually 

makes them interfere with each other. Resulting, the true evolution leads to an oscillatory behavior (in cosine) while 
this fictitious "classical" evolution gives a simple exponential damping (since dp^^/dt = A — 2pj^j^ and so on). 

Moreover, this evolution can be derived from a "classical" Poisson bracket, defined from our Moyal bracket {•, - jj, 
by dropping the Ty-translation from equation ([ll 



{/,5}d = ^ (20) 
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or equivalently 

{{f,g}ei)a = -^tyD^Y[F,G\Y. 

Due to these extra y-opcrators, one can check that this "classical" Poisson bracket does not satisfy the Jacobi identity 
unlike the Poisson bracket defined as the commutator of the 7k-product. Therefore, at the discrete level, we do not 
have a choice: only the Moyal bracket and not the "classical" Poisson bracket defines a proper Poisson structure 
(satisfying the Jacobi identity) and thus a consistent symplectic structure. 



IV. GOING BEYOND THE QUBIT 



Following Wootters' construction for discrete phase space, we consider a (/-dimensional Hilbert space H, with d>3 
a (odd) prime number (it can be generalized to the power of a prime number by working with finite fields), and we 
construct a.dxd discrete phase space as in the qubit case by introducing a suitable basis for the (i^-dimensional space 
of (Hermitian) operators on H. Choosing a basis of H, with vectors |o) labeled by a G Kd, we introduce the following 
two unitary operators: 

X\a) = \a + 1), Z\a) = Lo'^la), (21) 
where ui = exp(2i7r/rf) is a rf-root of unity. Then we define the dx d grid of operators: 

£»(p,g) = (22) 

a,o 

where the (finite) sum is implicitly taken over a,b € TL^. We remind the reader that 2 is invertible in TL^ and that the 

division by 2 is equivalent to the multiplication by the integer ((i+ l)/2. Let us underline the fact that this definition 
does not work in the qubit case for d = 2. Using the commutation relation Z^X"" = oj°'^X°'Z^ and the fact that co is 
a root of unity, we easily check that these operators are Hermitian and we can compute their traces: 

Dip,q)^ ^ D{p,q), irD{p,q) = l, tv D{p, q)D{r, s) = dSp^r^s- (23) 

We also check that J2p,q D{p, q) = dl. For any operator F on TL, we associate an observable / on the discrete phase 
space: 

f{p,q) = -^iYFD{p,q), F = f{P,Q)D{p,q), trF = ^/(p,g). (24) 

The constant function / = 1 corresponds to the trivial operator F = dl proportional to the identity. The operator 
F = Z gives f{p,q) = z{p,q) = ju'^ defines the space coordinate while the operator X'^ = X~^ with x{p,q) = jlo^ 
defines the momentum coordinate. The fact that the coordinates are given effectively by {io'^,uj^) instead of simply 
(p, q) hints that q- numbers for q=a; are likely to appear naturally in this framework . 
We define the ★-product as the operator product: 

{f^9){p,q) = ^trFGD{p,q). (25) 



q-numbers are usually introduced when dealing with q-deformed quantum groups. They are defined as: 




which converges to n in the classical limit q —* 1. 
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This ★-product is associative and non-commutative and its unity is the constant function / = 1. We define the 
corresponding Moyal bracket: 

{/, ?) = - ^ (/ * 5 - 5 * /) = - ^ tr D{p, q) [F, G\ , (26) 

where the Planck constant h is inserted by hand. We can give a more exphcit expression for the ★-product referring 
solely to the functions on space phase and not to the operators: 

{f*9){p,q) = ^ E fir,s)g{t,u)tyD{p,q)Dir,s)Dit,u) = ^ E f{P + r,q + s}g{p + t,q + u)co^^''-'-"\ (27) 

r,s,t,u (r,s),(t,u) 

which we compute by expanding the D-operators in terms of the basis X"Z^. We recognize the ★-product for the 

non-commutative torus with deformation parameter 9 = d/2Tr. 
This provides to a simple expression for the Moyal bracket: 

{f,9}*{P,(l) = ^ f{p + r,q + s)g{p + t,q + u)sm^{st-ru). (28) 

We can apply this formulate to the space and momentum coordinates, z{p, q) = ^w'' and .t(p, q) = and compute 
the canonical bracket (which actually corresponds simply to the commutator of Z and X): 

2 TT 

{z,x}^ = --ZX sm-. (29) 

Setting = 2 in this formula, we have sin7r/2 = 1 and we recover (up to a constant normalization factor) the Moyal 
bracket of the qubit case given earlier. We can compute the bracket of more generic functions: 

}^ = -| sin ^(a,5 - (.hi). (30) 

We would like to compare this to the standard Poisson bracket of classical mechanics in the continuum, {q, pj^t = 1. 
More precisely, we consider holonomy-like observables ^exp(iAq) and cxp(iAp). A simple calculation yields: 

{p»Aq^ e*^P},t = -A^ e*^i e'^P, (31) 

which is very similar to the formula above computed in our discrete setting. We can push this comparison further 
and compute the Moyal bracket of these holonomic observables. The standard Moyal product for a one-dimensional 
system is given as a power series in /i by : 

= (|)"E(-i)' (^) (^X"'/)(5r'5^5), (32) 

and the Moyal bracket is given by the commutator: 

{f,9}*st = -^{f *stg - g*st f) — > {f,g}st- 
in ft— »o 



Then the Moyal bracket of the holonomic observables gives: 



{e'^^, e'^^Ut = -\ e'^'i e'^^ sm^{a6 - (3j). (33) 



For dimensional purposes, since space and momentum coordinates do not iiave the same physical dimension, it would be better to 
consider observables exp(iAiq) and exp(iA2p) with the constants Ai and A2 a priori different and independent. In this case, all the 
formulas given here are still valid with A^ = A1A2. 
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Actually this is exactly the same as our discrete Moyal product if the following relation between the dimension d and 
the coefficient A is assumed: 

This is equivalent to lu = exp(2i7r/d) = cxp(i^A). This means that the discrete phase space structure is simply a 
truncation of the standard continuum phase space to specific holonomic observables defined above. Then the discrete 
Moyal bracket defined through the (finite dimensional) matrix commutator is exactly equal to the continuum Moyal 
bracket. Let us emphasize then that we have correctly considered the matrix commutator as defining the Moyal 
bracket and not an approximated discrete Poisson bracket. 

We have understood the relation between the discrete Moyal bracket and the continuum Moyal bracket, thus 
providing a representation of our discrete Moyal bracket in term of the continuum differential calculus. Instead, we 
would like to be able to express the discrete bracket in term of a discrete differential calculus, just like in the qubit 
case. To this purpose, we look at the translations. It is very easy to check X generates shifts in q while Z leads to 
shifts in p. More precisely, we have: 

X'Z^ D{p, q) {X'Z"-)^ =D{p + r,q + s). (35) 
Thus we can define finite difference operators as in the qubit case ^: 

idpf){p, q) =^ tvD{p, q){Z^FZ -F) = ^ tvF{ZD{p, q)Z^ - D{p, q)) = f{p + 1, g) - /(p, g), (36) 



and 



{5,f){p, q)^\ti Dip, q){X^FX -F) = \ir F{XD{p, q)X^ - D{p, q)) = /(p, q + 1) - f{p, q). (37) 



We would like to express the discrete Moyal product/bracket in terms of these differential operators and compare 
it with the standard formula of the continuum limit. First starting from equation (j27p . replacing lj — exp(27r/(i) and 
expanding the exponential, we find a formula similar to the usual formula (j32p : 



(/ ^ 5)(P, 9) = ^ E ^ (^) " E(-l)'^ f " E -'^"^V(P + r, g + .) E ^"-'^'^(p + t,q + u) 



k=Q 



This suggests that l/d plays the role of an effective Planck constant, as already implicit in the formula (I34p above, 
but does not allow to express the product in term of the differential operators. On the other hand, we can compute 
explicitly the iterated action of the differential operators: 

{S',S'j)ip, g) = E E(-l)'-''(-l)'-^ ^ ^ ^ j /(p + r, g + 5), (38) 
which is easily inverted for r, s > 0: 

/(p + r, g + .) = E E ( ) ( ' ) iW)(P' (39) 



We could also define differential operators of the other type such as (dpf){p, g) = ^ tr D{p, q)[Z, F] and the same for dq in term of X. 
Their action is not as simple as the other operators Sp and Sq and seem to involve shifts in p and g by ^ , which is harder to interpret 
on the lattice. Indeed, keeping in mind that z(p,q) = 2'^'^' ™^ compute 

<?) = ^ D{p, q)[Z, F] = (a;" * / - / ^ = -c^-J (^/(p + 1, g) - /(p - 1, g) 

Let us insist that ^ must be understood as an element of and equal effectively to the integer [d + l)/2. Indeed P + | does not 
correspond to a point in the lattice close to the original point (p, q) but to a point on the opposite side of the lattice. 
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This ahows an expression of the Moyal product in term of the (5-differentials, but we have not been able to simplify 
further in order to make it look like the standard formula. 

We now move on to the dynamics on the discrete phase space. Let us on the simplest dynamics: the "free motion" 
generated by the Hamiltonian H = {X + X'^)/2 which is the equivalent of the standard . More precisely, we 
can't take H = X since the operator X is unitary but not Hermitian. The eigenvalues ^'^of X are lo^ while the the 
eigenvalues of H are cos . Moreover we can compute the Hamiltonian function on the phase space lattice: 

Thus, besides the pre-factor ^ and a constant shift, this choice of Hamiltonian gives the standard p^ dynamics in the 
large dimension limit c? — > cx). The Hamiltonian evolution of the density matrix dp/dt = —j^[H^p\ is easily translated 
in term of the probabilities p{p, q) on the phase space using the Moyal bracket : 

Let us remind the reader that i is considered as an element of TLd and equal effectively to the integer {d+l)/2. Indeed 
g + ^ does not correspond to a point in the lattice close to the original point (p, q) but to a point at the opposite side 
of the lattice. We also point out the pre-factor which depends non-trivially on p. 

We can check that the probabilities q) for fixed p are conserved. These correspond to the classical trajectories 

-motion along the g-axis with fixed momentum p- or equivalently to the eigenvalues of the Hamiltonian H . This could 
be seen directly from the definition of the /^-matrices. Indeed, for fixed p, the operator 

^ ^i?(p,g)-^^f"X'^ (42) 

q a 

is the projector on the eigenvalue uj~^ of the (unitary) operator X. It commutes with X and thus with H . Therefore, 
we have: 



E 



= -\trV,[H, p\='- tr[H, V,]p = 0. 



Finally, one can check that such an evolution leads to negative probabilities (while keeping of course the sum of all 
probabilities normalized to 1) even if the initial probabilities were all positive. 

We could introduce a concept of classical evolution by using a simple discrete symplectic bracket: 

{f,9}ci{p,q) = -^[^pMs - dqfdpg], 

in order to define dynamics which would keep the probabilities positive under time evolution. This could be considered 
as a first order approximation of the full discrete Moyal bracket. However such a bracket does not satisfy the Jacobi 
identity and thus does not define a proper Poisson structure. At the end of the day, the discrete Moyal bracket {•, - j^, 
seems to be the only consistent choice for a discrete Poisson structure. 

One can find a rigorous discussion of "classical" vs quantum states in [1, [1] . "Classical" states are defined as having 
a positive probability distribution (for all choices of Wigner functions based on a fixed set of mutually unbiaised basis). 
Further physical explanations of why such states can be called classical in term of computational speed-ups can be 
found in Unitaries that preserves this positivity criteria, thus sending "classical" states onto "classical" states, are 
shown to form a subgroup of the Clifford group. The Clifford group being discrete, any continuous unitary fiow would 
necessarily produce negative probabilities and non- "classical" states. Nevertheless, for a fixed (time-independent) 



The eigenvectors of X are [a) while the eigenvectors oi H = {X + X^)/2 are sin{2npa / d) \a). 
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Hamiltonian, if we call r the first time at which the evolution is given by a Clifford group element (if it exists), then 
at all times nr with n G N, the evolution would again produce a Clifford group element. Thus, it could be that the 
evolution looks classical if we look only at the system in term of a fixed time unit. Nevertheless, this last remark is 
vague and remains purely speculative. 

Finally, we look at canonical transformation, i.e changes in the space and momentum coordinates that preserve the 
symplectic bracket. Up to now, we have been working with z = ui'^ and x = ui^ - we have dropped the ^ factor in front 
of the coordinates since it does not affect the bracket (P^. Let us define the new coordinates, with a, /3,j,S G 'Zd- 

i = ^"9+/3p^ J ^ ^79+<5p^ oj. equivalently: 
As we have already computed previously, we have : 

2 TT 

{z, x}i, = ——zx sin -^{o^S — f3j), (43) 

which gives the same bracket as initially as soon as {aS — f3'j) — 1, i.e the determinant of the change of variable 
{p,q) — > {p,<j) is 1. This is the behavior expected from classical mechanics: it ensures that the map {p,q) — > ip,q) 
and thus {z,x) — > {z,x) is one-to-one. As shown in j^], such a canonical change of variable actually corresponds to a 
unitary change of basis up to a phase, i.e there exists a unitary matrix U depending on a, (3, 7, 5 such that: 

UD{p,q)U'^ = Lui^P'^-P'^^ D{p,q). 
This concludes our analysis of the symplectic structure of these discrete phase spaces. 



Conclusion 



We looked at Wootters' construction of discrete phase spaces and discrete Wigner functions for quantum mechanics 
on finite-dimensional Hilbert spaces. We studied the induced Moyal product and the resulting discrete symplectic 
structure. Considering these phase spaces as simple examples of non-commutative geometries, we defined a discrete 
differential calculus and expressed the Moyal bracket in terms of these differential operators. We showed that the 
Moyal bracket satisfies all the same properties as the usual Moyal product in the continuum. Actually, we further 
proved that, for odd prime dimensions, the continuous Moyal bracket evaluated on certain holonomic observables 
gives exactly the discrete Moyal bracket. Resulting, the discrete phase space appears to be an exact truncation of the 
standard phase space. 

We also discussed the dynamics induced by the discrete Moyal bracket. The Hamiltonian evolution naturally leads 
to negative (pseudo-) probabilities, which is the signature of the quantum regime. This is expected since the Moyal 
bracket was constructed to describe exactly the quantum evolution on the finite Hilbert space. We show that one could 
construct using the discrete differential calculus a kind of classical Poisson bracket for which probabilities would stay 
positive. However, it turns out that such classical Poisson bracket does not satisfy the Jacobi identity (associativity) 
and thus does not define a consistent symplectic structure. This means that the discrete Moyal bracket remains the 
only consistent symplectic structure at the discrete level: we are in the deep quantum regime with no equivalent of 
the h ^ limit and no analog of the Poisson bracket. 

Finally, we only considered discrete two-dimensional phase spaces. It would be interesting to generalize our analysis 
to higher dimensional cases with a quantum system evolving in more than one space dimension. In fact, ,we can simply 
take the tensor product of the discrete phase spaces considered here. This is actually the Wootters' prescription for 
Hilbert spaces with non-prime dimension: if the Hilbert space dimension d is non-prime, we consider its prime 
number factorization d = and we construct the phase space as the tensor product of the discrete phase 

spaces associated to each prime factor d"' . This way, we get a 2n-dimensional discrete phase space. Nevertheless, it 
would be interesting to obtain a construction where the space dimensions are more intertwined with each other such 
that the space itself is endowed with a non-commutative structure. Maybe, this is given by the applying the same 
construction for non-prime dimension d than for a prime dimension without using the tensor product of its prime 
factors. We would then get a natural working example of quantum mechanics on a finite non-commutative space with 
a consistent symplectic structure and Moyal product. 
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